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written in APL are proposed for studying ways of dealing 
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The effectiveness of many systems depends upon the 
Begree Ot the calibration of their subsystems. For example, 
a ship with navigational equipment that is out of calibration 
may not be able to locate its destination, or, in the case 
of a Navy ship, locate an adversary. If the Navy ship's 
feapen syseem 15 also out of calibration the difficulties 
are compounded. An analogous problem arises in connection 
with engine de-tuning, when fuel consumption will likely 
increase and performance decrease, and with drift of communi- 
Cation systems. The detrimental effect of mis-calibration 
is well recognized: Navy ships and other systems are taken 
meebanges or Other facilities for testing and re-calibration. 
The purpose of this thesis is to set up mathematical 
models for studying ways of dealing with mis-calibration. If 
the various aspects of the problem can be assembled, some 
mecance 15 then available for dealing with it effectively. 
Although various realistic elements of the problem can be 
introduced, the fundamental issue is this: given that 
important subsystems depart from calibration and effective- 
ness as time passes, it is desirable to determine a schedule 
men re—-Calibration that (nearly) optimizes system operational 
effectiveness. Frequent calibration of important systems 


would be highly desirable if this were a cost-free operation, 


but in reality the operational cost of calibracwonesasea ee 
time during which the system is unavailable for, or so 
degraded as to be incapable of adequately performing, its 
Operational purpose. Figure 1 is an idealized graph cE 
operational effectiveness against time. The periods of 
duration C denote those periods during which the system has 
zero effectiveness because it is undergoing calibration and 
hence is out of the operational area; the periods of duration 
T represent those periods during which the system is opera- 


tional, but of diminishing effectiveness. 


e(t) 





ee CL 
a C T 


©) 


Figure 1.1. Idealized Graph of Operational EEfGeCtiy ei 


a 


The graph suggests that if effectiveness drops with time 


Pere will be an Optimal value for T, a "best" period, T*, 


at which to calibrate. We now show how such a period may 


be determined. Later, more complex and realistic models 


and simUlation results will be introduced. 


de: 


Il. MATHEMATICAL MODELS 


For a mathematical treatment let e(t) be the effective- 
ness, e.g., the probability of successful mission completion, 
at time t after the calibrated system returns €6 Senpviecs 
Let C be the time required for calibration, and T the duty 
Or on-station time. Then the average effectiveness over a 


cycle of length T+C, and hence in the long run, is 


f e(t)dt +0 
e(T) = SSC; (22) 


wae 
the term O represents and emphasizes the total lack of 
effectiveness during the calibration period. In order to 


maximize e(T) it is useful to study the derivative 


a 
(T+C)e(T) - f e(t)dt 


dee) 0 
SSN ef ee (22525 
aT (p + cy 2 


as it depends on T: if de(T)/dT = 0 for T* > 0 then sta 
a candidate for a time between the end of one calibration 
and the beginning of the next. Equivalently, (2.2) asks if 


there iS a positive solution T*, of 


a 


e(T) ara of elt)dt (2.3) 
0 


eZ 


for fixed positive C. The fact that such a solution always 
exists, and that it defines an optimum can be established 
from the usual second derivative Sree aA Since the 
Setimal £ satisfies (2.3), it turns out that at the optimum 
the average effectiveness over an entire cycle equals the 
effectiveness at the time the active part of the cycle 


ends; or symbolically 
e(T*) = e(T*) (ody) 


where the over-bar signifies the time average of effective- 


ness over T* +C. 


To build understanding, examine some extremely simple 


specific models. 


eo INEAR EFFECTIVENESS LOSS 


Eel te 


Ze) 


Semetat the downward-sloping parts of the graph of Figure 
imeeeare strictly linear. Then (2.3), the equation for 


optimal T = T*, is 


eer = — (T-5T), O0< T<a ; (2.6) 


gs 


1 


It is clear that no value of T > a ~ can be Gptamunee 


equation (2.6) simplifies Go. emer euadrauus 
2 
aT* + 2aCl —.2 Ga =a0 ( 2.9— 
with a single positive solution 


T*¥ = -C + Vc? + 2C/a (2288) 


at which the optimum value of effectiveness 





Tx 
J e(t)dt 
Re €- at) Z iS ee pas cf eee 
Se) a= eae eae perce = 1 + aC. \ a5Gae a6 
= ese a eee (2.9) 


It is interesting that the solution depends only upon the 
parameter aC the product of calibration drift rate, a, and 
the length of the re-calibration period, C. For instance, 
if ac > 0 then effectiveness approaches unity if either the 
rate of calibration degradation, a, approaches zero, oFwea 
calibration time, C, approaches Zero) onus bormy cine 
approaches zero more rapidly than the other gets large. 
Alternatively, this shows that equal-effectiveness or a-C 


tradeoff curves are simple hyperbolas in the (a,C) plane. 
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The above model is rather crude, but is easy to under- 
stand. There follows another model that 1s more gqualita- 


Eively appealing. 


B. LINEAR DEGRADATION WITH DIFFUSE DAMAGE 

Consider next a more specific model for effectiveness, 
one that relates to damage inflicted on a target after time 
t has elapsed, and the system has developed an (unsuspected) 
bias of magnitude at. At that time the x-y error made in 
locating a target is assumed to be given by the joint 


Gauss/normal density 





2 
eee = =~ exp ea 
G 


If a cookie-cutter damage function with radius R is in 


effect (no damage if x? + y* > Re destruction if 
x? re a Peo cen 
NG a coc, een y 


(x* + y? < RY) 
However, this is difficult to work with, and even overly 
Simplistic. Instead, suppose that a von Neumann-Gauss 
diffuse damage function can be used; i.e., that the proba- 
babedty Of Critical damage to a target located at (0,0) by a 
weapon with impact point (x,y) is equal to 6(x,y) = exp(-a(x 


Me-)). Then 


eS 


av) 
ain 
II 


[f 6€x,y) £(x,y7t) dxdy 








= a f i exp(-a (x+y) ] £(x,y;#) axdy 
2G 9 =o —< 
ce 2 2 
= ( - : exp(- 58°78) exp [-ax]dx) (2.1 
V2NG =o 0 


by virtue of the symmetry assumed; almost free of charge we 
can consider asymmetrical damage functions, but the oppor- 
tunity is declined. The above integral is evaluated at sight: 


it is seen to be essentially the convolution of two normal 


densities. After squaring, as demanded by (2.11), 
lie Gee 
e(t) = —- + exp [- re (2 ee 
(o° +1/2a) (o~ +1/2a) 


Instead of dropping off linearly, as in the previous Gase, 
e(t) first diminishes rather slowly, later falling quite 
rapidly (exponentially fast) towards zero: by the time 


at = Vo* + l/f2a, effectiveness is just below 40% of its 





maximum, while if at = 0.5 yeee l/f2a, effectiveness i585 about 
78% of the maximum; finally if at = 0.25 ao + l/2a, eitecs 
tiveness is 94% of the maximum. Note that the maximum 
effectiveness is (l + Die le < djeat southies a or a become 
large, meaning that if either weapon effectiveness falls off 


rapidly with miss distance (a large) or the ultimate weapon 
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; : Z . 
Pewee mar iciicemis meat (o 9 large), then even maximum 
effectiveness is low. 


inorder to sSolvyemror the optimum T* write 


L 


(T+C) exp(-(aT)7/(074+1/2a)) = f exp(-(at)*/(o7+1/2a)) dat 
0 
eee 
Change the variables to the dimensionless version 
Ty cy) a eee = (aC) (c-0720)- , (2.14) 


so one can solve the following dimensionless equation once 


eae tor all for t*: 
2 Z 
(t+y)exp(-1t~) = f exp(-z")dz ; el) 
0 


the positive value of t, namely t*, that satisfies this 
equation may be located by Newton-Raphson, or even graphically: 


eme, can plot, for given jy, 
2 
imma se (tty exp (—7 —) 
and 


Cpa i exp (-z~) az 
0 


on the same piece of paper, vs. tT. 


iy 


The arbitrary selected y values and the corresponding 
T* values from the computer program which solves the dimen- 


Sionless eguation (2.15) are shown in Table l. 


TABLE Ll 


Gamma and Optimum Tau Values 


GAMMA TAU 
0.001 On Gigs 
02002 OS 
0.003 0.144 
0.004 0.165 
0.005 0.181 
0.006 0.195 
On007 0.207 
0.008 0.218 
0.009 OF 27 
0201 Deze 
(02 es 02 
O02 0.346 
0.04 0.381 
00> 0.410 
0.06 0.434 
0.07 0.456 
0.08 0.476 
0.09 0.494 
Oe O4oaMl 
0.2 fic s2 
Gis ene 
0.4 Oenime 
0-5 0.825 
0.6 0.867 
G7 0.904 
0.8 0.937 
0.9 0.966 
diye 0.992 
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In order to redefine t and y in (2.14) and the effective- 
ness formula (2.12) in more meaningful form, again change 
variables to 


ee ee ee ke = Caso (2.16) 


a +1/20 
where v and p might be called vulnerability and probability 
pmesuiccess respectively, and k is constant. In the simula- 
tion chapter connections are developed between v and the 


radius of a (roughly) equivalent cookie-cutter damage 


mumction. So one can write 
T =k vl-p T ; y¥ = k vl=pc. (eo a7) 


We focus attention on the representation (2.12) in what 
follows, mainly for analytical and computational 


convenience. 
e(t) = p exp i (eee C27 Ss) 


Thus, the preceding expression at optimum leads to the 


relationship 


e(T*) = p exp (ert) (2 219) 


muemconsequently, the optimal proportion of on-station time 


can be obtained as follows: 


i 


ej 
> 
rH 
+ 


* 
rt (2.200 


k VY1-=p 

so the proportion of on-station time 1s, Undervepeimum 
COndlLiLens, 

arte = ory (2. 20m 
Since optimum t values are available from Table l, one can 
very easily calculate the effectiveness given some constant 
variance and v or only p. Some of the results are tabulated 
in Table 2 as an example. 

Additionally, plots of tau vs. gamma and effectiveness 
vs. tau are presented in Appendix A. It is observed that the 
effectiveness decreases as the variance increases while v 
is held constant. Effectiveness vs. tau plots illustrate 
the behavior of the effectiveness representation in a more 
understandable fashion than does the formula itself. 

Example: Suppose a = 1.5 yds/month7.@ — 02 oe mene 
57 = 20 (yds)? and p = 0.9 are given. First find ygeeon 
(2.17) as 0.053, then look up corresponding t* value from 
Table 1 which is 0.417. Later, from (2.20) T* is 32793 epee 
and from (2.21) the proportion of on-Station time aseesaee 
and from (2.19) or Table l an average effectiveness of 


75.6% can be obtained. 


20 


TABLE 2 


bmepectiveness fOr ConstanteVariance and v 


EFFECTIVENESS EFFECTIVENESS 
(o* =10; v= 200) (o* =20; v =200) 
OS 2 vA20s 
0.940 Oreo 7 
Oe Ome2 0 
OT Necc> 
Oe g22 0.880 
Or oelay, Oro.2 
Deg i OG 7. 
0.908 Oe cio 
0.904 O63 
Or 9 Oi 0.860 
Ore 69 Ores 0 
0.845 0.806 
Os. 4 0.786 
G7 305 Oe oS 
Oecd 02753 
Ode On eC 
Oo O25 
0.746 Orb 2 
0 .- 7238 O00 
0.639 Oe Eu 
Oro. OF S47 
Orns 272 0.499 
0.482 0.460 
0.449 Orage 
0.421 0.401 
OS 26 0237 3 
0.374 Oe 
356 0. 340 


Zt 


C. LINEAR DEGRADATION WITH DIFFUSE DAMAGE USING RANDOM DRIFT 
For an alternative model, that incorporates the woes c esas 

different drift rates of different individual ships or 

systems, suppose that the drift, a, 1s a random variable 

with an appropriate distribution function inseceddmeme 


Constance as im (Zee namely, the effectiveness conditional 


on a is 
Z 
(aie) 
cc.) = ee (2.2m 
. (o° +1/2a) (o” +17 2) 


Then, the expected average or unconditional effectiveness 


over a cycle of length T+C, in the long nun? eeccomes 


lt 
f Ele(t;a)]dt 


— eet 
Ele(T;a) J 7 Saar (2.255 
in order to be specific (but not necessarily realistic) and 
also so that explicit mathematical results are obtained, let 


ao have a gamma distribution function with parameters A and 


By 
—-AX B-1l 
fa-(x;A,8) =) A2 ee ee (2.205 
. Ee) 
and put for fixed ae i.e., the square of the deifrt mae 


away from Calibraugren, 
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4 
Ee) = at = -* so were a ol (202 OD) 
= A fe 12 B 
mew Use (2.22) and (2.23) to obtain 
Woe * ue 2 
Ele(t;a)] = 5 Ole eee) | bae (074 ,6) dx 
- (omer 17 20) 80 G6 +tl/2a : 


after substituting the gamma density function it is easily 


seen that the result of the integration yields 


2 


Q 
Ble(t;a)) = —,—@* {A ___ (2.26) 
2 (o° +1/2a) ae fz 


anyon 


Or equivalently, in view of (2.25) 


B 
Bie (t:a) | = —— 5. (C72 7.) 
. (o° +1/2a) js (at) 
2 
(o +1/2a)8 


Various analytical properties of the previously described 
model will now be recorded. These provide useful insights 
into the behavior of the effectiveness at time t. 

mee if 8B + wo, (2.27) becomes 


ee 2 


Beer; at) | 5 Selle (ac) C2526) 
o0 +1/2a 
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which reflects the fact that if 8 increases the variance of 
a in the distribution of drift rate decreases towards zero, 
and the situation reduces to Ehat of Modelee 


Zee Le be lee Enen 


de 
2 ‘ 


EB LesGe a) ia 5 
(o° +1/2a + (at) ~) 


(2.225 
which is larger than the effectiveness in the equal-drift 
case. 

In order to solve for the optimum T* for the general 


case write 


T+C 2 1 
(at)? BO 9! (at)? 5 
(1 +—, =" —_) (1 +—3*+__) 
(o +1/2a)86 (O° Wey 2a 
Change the variables to 
t = (aT)/(o74+1/20)17* ; y = (aC)/(o~ 41/20) een 
sO one can solve the dimensionless equation 
T 
iis | eee (2.328 
4t8 0 14248 
B 8 


the positive value of T, namely T*, that satisfies this 


equation for any constant 8 may be found by janecemnouees 
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program. Im fact, one may get the solution for the special 
case 8 = 1 by making use of arctg integration for the right- 


hand side: 
2 4 
We acecea ty lstoy ) = 7 « (Zo?) 


In general, the right-hand integral can be transformed 
to the integral of a Student's t density, and the t-tables 
found in most statistics books can be used to evaluate it. 

Again, the arbitrary selected y values and the correspond- 
ing t* values for g = 1 are presented in Table 3. 

At this point, it is very easy to calculate the effec- 
tiveness given some constant variance and v or only p from 
(2.19). Some of the results are listed in Table 4 as an 
example. 

In addition to the tables, plots of tau vs. gamma and 
effectiveness vs. tau are presented in Appendix B. As in 
the previous case, effectiveness decreases as the variance 


jfmaereases when v is held constant. 
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TABLE 3 


Gamma and Optimum Tau Values (8 = 1) 
GAMMA TAU 
OO Or 02025 
002 021955 
O.008 0.145 
0.004 OGG 
O00. O62 
0.006 0.196 
0.007 0-209 
0.008 0.220 
O00 Oie230 
0.02 e233 
OFe 2 0 307, 
OF 0s 0-355 
0.04 Oro 
O05 0.424 
0.06 0,451 
0e.07 0.476 
OOS 0.499 
0. 09 OF 520 
OF 1. 0.539 
eZ On6o7 
O.3 0.793 
0.4 0.879 
Hiss 0, oa 
0226 i Oaks 
Oi..7 LO 
0.8 tel 
On io 
a0 lae2s 
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TABLE 4 


Effectiveness for Constant Variance and v (8 = l) 
PREEECTIVENESS BRE ECT EVENESS 
(o* = 20; v =150) (o* =30; v=150) 

Ome 82 e338 
Deo iult Oi oe 
0.864 OeoLG 
nO a Orreaint 
0.854 0.807 
0.850 Ores 02 
0.845 eos 
0.842 vo 
0.838 02791 
O85 5 O77 38 
0.806 Om ol 
0.784 0.740 
OG O27 22 
0.748 0.706 
Os 5 ee 692 
ees 02679 
O27 06 0.667 
0.694 (ios 6 
0.684 0.646 
OR. 29 vn 506 
0.542 eo 2 
0.498 0.470 
On igs 0.437 
0433 0.409 
0.408 ies co 
Omiser 0. 3 oG 
Orso 0.348 
OR a2 332 


Ze 


III. TRANSFORMATIONS AND SIMULATION 


A. TRANSFORMATIONS 

Earlier it has been shown that tT* can be computed in 
terms of y. In order to simplify this step, 1€ would be 
desirable to be able to represent t* by some simple formula 
in terms of y. Following the lead of statistical regression 
studies, it iS sometimes possible to investigate the effects 
produced by transformations of the predictor variables; GE 
by transformations of the response variable, or by both. 
Clearly there are many possible transformations of gamma and 
tau values. Several different transformations of gamma and 
tau could be tried for the same model, of course. The choice 
of which is sometimes difficult to decide and the choice 
would often be made on the basis of previous knowledge of 
the gamma and tau under study. The purpose of making trans- 
formations of this type is to be able to use a simple 
regression model in the transformed tau and gamma, rather 
than a more complicated one in the original gamma and tau. 
Some suitable transformations of gamma or tau can also be 
found by plotting them in various ways. First, In(taujeigee 
In(gamma) has been plotted for the linear degradation with 
diffuse damage case and the linear degradation with diffuse 
damage using random drift case, later various power trans- 


formations have been applied to tau values and simple 
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regression equations have been derived in order to obtain 
the optimum tau values directly for arbitrary selected gamma 
values without having to go to tables or equations (2.15) 
and (2.32). Some of the transformation plots are shown in 
Appendix C. After obtaining a suitable power transformation 
of tau, one could guidely calculate the effectiveness values 
by using predicted tau values from Oe ion equation. 

This attempt at simplification deserves more study before 


it can be said to be truly satisfactory. 


Be oLMULATION 

Simulation is essentially a controlled statistical 
sampling technique (experiment) which is used, in conjunction 
with a model, to obtain approximate answers for complex 
(probabilistic) problems when analytical and numerical 
techniques are too expensive, or infeasible. 

The main purpose of the simulation in this thesis is to 
be able to evaluate effectiveness for other kinds of damage 
Pmimceions OF error distributions which are difficult to work 
with, as alternatives to a von Neumann-Gauss diffuse damage 
function. For example, if a cookie-cutter damage function 
with radius R is in effect then a closed form solution of 
the effectiveness similar to (2.12) is not as simple. Unlike 
a mathematical solution, the answer one obtains from a 
Simulation is an estimate of the effectiveness. It is 
absolutely necessary to have some idea of the precision of 


the effectiveness. For this reason, the effectiveness 
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estimated from each simulation has error bounds of two standard 
deviations for valid comparisons. The interactive simulation 
program, written in APL, 1S presented in Appendix D. 

The scenario developed in this program determines the 
optimal time for a submarine to come in to port for instru- 
ment re-calibration. In other words, simulation is being 
used to determine T* for a variety of cases--ones in which 
the previous neat mathematical theory of Chapter II cannot 
easily be extended. A vector of possible times (in arbi= 
trary time units such as days) at which the submarine should 
be brought back for equipment re-calibration is needed. For 
each of these times the program estimates the expected 
effectiveness of the submarine. The time that corresponds 
to maximal effectiveness is considered optimal. Although the 
effectiveness of the submarine changes continuously with 
time, in the simulation the effectiveness is estimated only 
at discrete but closely-spaced time points. The more perme. 
one has, the smoother the effectiveness curve, but the longer 
the program takes to run. In addition, the duration of re- 
calibration of the equipment (C) and the number of replica- 
tions of the simulation should be entered. Again, the 
precision of the estimates of the effectiveness curve gets 
better with more replications, but, again, it takes longer 
to run the program. 

Effectiveness is measured as the probability of damaging 


a target ship that is 1000 distance unlts away et reneeae 
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Submarine. The weapon is a straight-run classical torpedo 
with a proximity fuse. The submarine fires the torpedo along 
some bearing and the torpedo is supposed to explode at the 
Beihiitenmearest touthe target. But the equipment to locate the 
target develops calibration problems with time, namely 
calibration drifts by a certain distance for every time unit 
according to the following specific alternative model 
options, any one of which may be considered by the analyst: 

l. It might get deterministically worse with time. So, 
the expectation of drift becomes E(at) = at; the rate a must 
be specified. This is Model l. 

2. It might get randomly worse with time. On day T, the 
drift is mismeasured by an amount T x Normal(0,Sigma). 


Although the mean error is zero, the variance of the error 


increases as ence deine Notice that the random multiplier 
is constant in each replication of the simulation. MThis is 
Model 2. 


Peeereomight fluctuate randomly with time. On each day, 
mae drift 15 mismeasured by an additional drift error. This 
error term is random and comes from a Normal(0,Sigma) dis- 
tribution, where sigma is expressed in distance units and 
represents the standard deviation of the error distribution. 


Notice that the expected error is always zero, although the 


variance grows proportionally with time. In this case, 
calibration can improve or worsen with time. This is 
Model 3. 


Sek 


4. It might exhibit a random drift, with magnitude drawn 
from a gamma distribution. Thus, calibration drift forse 7m 
time unit becomes a gamma random variable with shape param- 
eters lambda and beta; the drift gets worse with time. Note 
that the program uses GAMMACK APL library function to 
generate the incomplete gamma random variable. This is Model 
4. The gamma variability explains the differences in drift 
exhibited by different system copies. 

The user may choose which of the above models best des- 
cribes his or her situation. 

After a straight-run classical torpedo is aimed at a 
point influenced by one of the preceding calibration errors, 
it may not explode at precisely the point on that bearing 
that is closest to the target, i.e., the proximity fuse is 
assumed to be not perfectly accurate, as is true in reality. 
The error between the closest point and the explosion point 
can come from either normal distribution with variances in 
the X and Y direction or a uniform distribution with (gee 
and -(=Y71)= 

Finally, the target is damaged with a probability calcu- 
lated according to one of the following optional funceren 

1. Gauss diffuse damage: Probability of critical damage 
to a target located at (0,0) by a weapon with impact point 
(xX,;Viets equal to>6 (x34 )=— exp (=a (x> yom 

2. Cookie-cutter: The torpedo will destroy the earcee 
if it 1s within a certain radius R, and it will do no cama 


1f it is outside this radius. 


a2 


Bee Trapezoidal: The trapezoidal damage function has a 
central, circular plateau of radius Ri - If the target is 


within R, distance units from the exploding torpedo, then 


1 
the target is damaged with probability one as it is in the 
cookie-cutter. In addition, the function has an outer 

eorcular rim, Of radius Ro, Ro > Ry, beyond which the proba- 
ali ty of damaging the ship is zero. Between the two 

radii, the damage probability goes down linearly. 

As a result, the simulation program provides three basic 
Sutput arrays. EPR containswthesestimated effectiveness at 
each time increment, delta, out to the maximum time. Again, 
the effectiveness is simply the probability of the torpedo 
destroying the target. STDEFF contains the standard devia- 
tions of the estimates in EFF. Lastly, AVGEFF contains the 
long-term average effectiveness of the submarine if it 
returns after the various times following calibration delays 
moult at the beginning of the program. 

Plots of simulation results in Appendix E were obtained 
from the general plot function in GRAFSTAT which is an APL 
workspace for the interactive creation of scientific-engineering 
graphics, for interactive data analysis and for the inter- 
active development of APL graphics output routines. It runs 
on both the IBM 3277GA graphics terminal and the 3278/79 
terminal. Full color control is available when running on 
the 3279 terminal. The simulation program is attached to 


GRAFSTAT in order not to waste time in assessment of the 
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results and to be able to make sensitivity analyses easily. 
For each run, two different graphs are obtained. One of them 
shows estimated effectiveness of a submarine with error 
bounds of two standard deviations. The other one demonstrates 
the long-term average effectiveness of a submarine, given 
that it stays at sea for a tour of length T. Siieestnow 
is a large number of combinations in the program, only 
limited numbers of results are presented in this thesis. 
The first six plots simply exhibit the sensitivity of the 
system effectiveness to the parameter alpha of the Gauss- 
diffuse damage function; other variables are held constant. 
It is observed that the optimum T gets smaller as the 
parameter alpha gets bigger; since a large a corresponds to 
a small effective damage region, a short tour length is 
necessary to keep effectiveness high. In addition, some 
interesting combinations of alternatives are also presented 
in order to give an idea about the behavior of the other 
parameters. There is not much to say about them since they 
are quite self-explanatory. 

It is of interest to check out the possible relationship 
between Gauss diffuse and cookie-cutter damage functions. 
An analytical-numerical solution for the optimum interval, 
T*, is available for the von Neumann-Gauss function, but 
none is for the cookie-cutter or the trapezoidal functions. 
If the latter can be reasonably matched to the former, an 


approximate analytical solution 1S avyatlables-onwaen-— 
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latter damage functions. Possible matches to cookie-cutter 
are these: 
00 R 
ieeicam matehing : f rexp(-ar~)dx ay ee) ees Rae 
So, the Gauss diffuse coated with parameter a ia cookie- 
cutter with radius R are approximately equivalent. 

Zee Median matchang: exp (-ar“) = 0.5 which is equivalent 
to a(R) = tee Wee Thus, this transformation matches the 
medians of the two functions. 

It is not possible to find a unique best transformation 
between a and R for all cases. However, one can test a 
proposed transformation for each specific case by simulation. 
For demonstration purposes, a(R) = In (5/3) /7R- 1S a proper 
Bieam@erOumation for R = 14.29 on condition that the other 
variables are held constant. It provides the same optimum 


T for Gauss diffuse and cookie-cutter damage functions as 


is seen in Appendix E. 
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IV 2 CONCEUS TOR 


The purpose of this thesis has been to show that mathe- 
matical models, augmented by a computer simulation program, 
can provide useful ways of studying the impact of miscali- 
bration upon operational effectiveness. We have concen- 
trated here on specific and convenient models, but it is 
obvious that other mathematical models can be treated 
Similarly. Other analyses similar to the ones we have 
discussed in the computer program can be conducted pertaining 
to the other alternatives. The relative effectiveness of 
different system configurations can also be investigated. 

The results of the simulation can then be analyzed to ascere 
tain the significance of different factors in various scenarios. 

The possibility exists that operational data will reveal 
different underlying distributions, and suggest alternatives 
for evaluating effectiveness other than the ones described 
in this thesis. The present thesis is to be considered a 


pilot study of the calibration issue. 


36 


PLOT OF TAU VS GAMMA 


PLCTS OF TAU VS. GAMMA AND EFFECTIVENESS VS. TAU 
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PEOLrS OF TAU V5. 


PLOT OF TAU VS GAMMA (A IS GAMMA DIST. WITH BETA=1) 
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LINE FITTING FOR LN(TAU) VS LN(GAMMA) 
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APPENDIX D 
COMPUTER SIMULATION PROGRAM 
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